Interacting Quintessence, the Coincidence Problem and Cosmic Acceleration 
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Faced by recent evidence for a flat universe dominated by dark energy, cosmologists grapple with 
deep cosmic enigmas such as the cosmological constant problem, extreme fine-tuning and the cosmic 
coincidence problem. The extent to which we observe the dimming of distant supernovae suggests 
that the cosmic acceleration is as least as severe as in cosmological constant models. Extrapolating 
this to our cosmic future implies terrifying visions of either a cold and empty universe or an explosive 
demise in a "Big Rip." We construct a class of dynamical scalar field models of dark energy and dark 
matter. Within this class we can explain why supernovae imply a cosmic equation of state w < — 1, 
address fine tuning issues, protect the universe from premature acceleration and predict a constant 
fraction of dark energy to dark matter in the future (thus solving the coincidence problem) , satisfy 
the dominant energy condition, and ensure that gravitationally bound objects remain so forever 
(avoid a Big Rip). This is achieved with a string theory inspired Lagrangian containing standard 
kinetic terms, exponential potentials and couplings, and parameters of order unity. 
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There is mounting evidence for the presence of an enig- 
matic dark energy component of the Universe. The na- 
ture of this component is one of the most profound mys- 
teries of physics (see 0] for a recent review). Dark en- 
ergy was first implied by a measurement of the luminosity 
distance-redshift relation, di,(z), through observations of 
supernova la (SNIa) at cosmological distances [2j,[3|. Re- 
cent SNIa data, galaxy redshift catalogues and analysis 
of the cosmic microwave background anisotropy favor a 
cosmic dark energy equation of state (wq) close to or 
even below -1 [111 110. 

Models in which this missing energy is a cosmological 
constant A fit this data well. However, these models have 
two serious drawbacks: the fine-tuning problem (A must 
be fine-tuned extremely precisely in the early universe 
relative to the energy scale at that time); and the coin- 
cidence problem (A must be set to an extremely small 
value in the early universe relative to the energy scale at 
that time such that matter domination and structure for- 
mation could occur) . Cosmological constant dark matter 
(ACDM) models are therefore not an explanation but at 
best an effective description of the dark energy. 

Simple models of dark energy as dynamical, non- 
interacting scalar fields, "quintessence" |8(, do not re- 
quire fine-tuning of energy densities or field values at 
early times, a because the field follows an attractor so- 
lution as it rolls down a potential [1E3- However, the 
potential must still yield a ratio of matter to quintessence 
energy density of pmo/pqq ~ 1/2 today, while allowing a 
matter dominated phase so that structure can form. Such 
models therefore do not address the coincidence prob- 
lem. Furthermore, there is difficulty achieving a current 
quintessence equation of state < —0.8 while retaining 
significant matter density [10J, except with quintessence 
potentials designed specifically for that purpose [Tl|. 

The possibility that wq < —1, i.e. that the sum of 
pressure and energy density of the dark energy is less 
than zero, has led to the idea of "phantom energy" [T^| . 
Explicit scalar field models of phantom energy produce 
super-exponential expansion of the universe by introduc- 



ing a negative kinetic term in the Lagrangian. The scale 
factor reaches infinity in finite time and the universe ends 
in an explosive "Big Rip." Setting aside theoretical dif- 
ficulties of such models, e.g. the fact that they can lead 
to unstable solutions by virtue of violating the Dominant 
Energy Condition [l3j . the cosmic coincidence is not ex- 
plained. In fact we would argue that the fine tuning prob- 
lem is exacerbated compared to cosmological constant 
models. Instead one appeals to an anthropic argument: 
if conscious life appears at all it must do so within the 
brief era of cosmic structure between the onset of mat- 
ter domination and the beginning of cosmic acceleration 
since the universe is destroyed soon after 0] . 

We look upon these difficulties with simple models as 
opportunities — guideposts in the vastness of theory-space 
to direct our attention. A natural expectation in the con- 
text of string- inspired cosmologies are interacting fields in 
the dark sector (see and references therein). Whilst 
the simplest models are tentatively ruled out by compar- 
ison to cosmic microwave background observations |16| . 
previous studies did not consider the possibility of more 
than one dark matter component. In this letter, we gen- 
eralize the class of interacting dark matter and dark en- 
ergy models suggested in |15l Il7| and demonstrate by 
means of an explicit example that this resolves the above 
difficulties. 

We call this class interacting quintessence cold dark 
matter (iQCDM) models. Our generalization consists 
in imagining additional field content of the dark sector. 
We propose that there is a quintessence-like rolling self- 
coupled scalar field <f> as well as other fields that can act as 
dark matter. These dark matter components will gener- 
ically be coupled to 4>, but the details of the couplings 
would depend on the details of the underlying funda- 
mental theory. In general one of the components will 
have a more <f> dependent coupling than the others during 
the brief period of cosmological history between matter 
radiation equality and today. We call this dark mat- 
ter component the interacting cold dark matter (iCDM) 
and group all the other components under the term non- 
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interacting dark matter (nCDM). 

Let the quintessence field <f> be contained in an effective 
potential that is the sum of two parts: a self-interaction 
term Vq(4>) = f (<j>), and interaction with the iCDM, 
Vm (4>) — g {4>) e~ 3N . The exponential decay (N is the 
logarithm of the scale factor) is a consequence of the di- 
lution of the dark matter with expansion, an important 
fact we will return to below. At early times the form 
of the potential should prevent quintessence domination 
through an attractor solution (tracking), while at late 
times <p will be at the minimum of the potential, which 
should yield a supernovae-measured wq < —0.9. This 
yields the following conditions: 
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where N is expansion e-folds and k = \/8ttG. These 
conditions arise from: a requirement that tracking solu- 
tions exist [HI, that the tracking solutions keep 4> sub- 
dominant [9j and that once at the potential minimum <j> 
will grow to dominate with an apparent wq < —0.9. This 
apparent wq is what supernovae measurements would fa- 
vor given a measurement of £Imq based on the assump- 
tion that matter did not interact with the quintessence. 
Steeply curved functions are favored by the second and 
third conditions. An additional criterion is that a signif- 
icant period (>7 e-folds) of matter domination precede 
the onset of quintessence domination. For a specific im- 
plementation one would like to find simple forms for / 
and g which satisfy these conditions and arise naturally 
in the low-energy limit of string theory. While this may 
seem to be a tall order we will now go on to show that 
we can simultaneously satisfy all these constraints using 
exponential potentials and couplings. 

As an illustrative example, consider the following spe- 
cific model: 
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where the quintessence field <j) interacts with a scalar 
iCDM field x with coupling constant 7, n = \/8ttG, 
C>nCDM contains nCDM, and Cb contains baryons. The 
constants a, j3, 7, tuq, 4> c would take values from an 
underlying physical theory. The value of <j) minus some 
offset cf> c sets the mass of x- The magnitude of the mass 
scale niQ is not actually important — it can be changed 
to current quintessence scales by a shift of <p and 4>c- 

One can "integrate out" x if it oscillates on a much 
faster timescale than the change of <fi. This calculation 
[22| results in an effective potential V, including the ex- 
ponentially decaying interaction term, and the following 
equation of motion for </>: 



_^3iTf + ^ = V^V Q + V M 

Vq = tyiqG qk (<W (quintessence) 

V M ee ^ g^ ^' (all matter) ' 

I ( M J + V Q+ V M ( total ener Sy density) 
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(Nq — N) is the number of expansion e-folds before the 
present. The effective Lagrangian leading to these equa- 
tions resembles that of exponential VAMP models |18j . 
except for the crucial difference of including additional 
species of dark matter (nCDM). We effectively model 
the impact of the baryonic matter and nCDM on the 
evolution of <f> through the term r/e^ K ^ 0_ ^ c ', which sets 
the ratio of this non-interacting matter to iCDM today. 
From the point of view of <fi, non-interacting species of 
matter have the same effect as the offset <j) c in the in- 
teraction potential. Thus, without loss of generality, one 
can absorb the magnitude of r into a shift of <p c . For 
simplicity we shall set r = 1. 

A schematic representation of the evolution of <fi in its 
effective potential V is shown in Fig. ^ The potential 
V is the sum of a positive and negative steep exponen- 
tial. Between them is a non-stationary minimum. We 
show the evolution of the fractional densities of radia- 
tion, matter, and dark energy for a numerical example 
in Fig. [21 At early times, <f> will be on one side of the 
potential, high above the minimum. To the left of the 
minimum, cf) will track down the potential, always re- 
maining a fixed, small portion of the total energy density, 
p<f>l ptot — 3 (wbg + 1) /a 2 0j until it reaches the mini- 
mum. The background equation of state wbg is \ dur- 
ing radiation domination, and after matter dominates. 
Primordial light-element abundance measurements con- 
strain the total energy density at the epoch of Big Bang 
Nucleosynthesis, and thus we must require a > 5 [l9j. 
On the right-hand side of the minimum, <fr will similarly 
track down this decaying potential to the minimum, keep- 
ing a fixed, small portion of the total energy density: 
P<t>l Ptot = 3wbg/20 2 H3 until matter domination or un- 
til it reaches the minimum. At early times this tracking 
keeps quintessence largely independent of initial condi- 
tions and protects the universe from premature acceler- 
ation. Note that iCDM domination is also not possible 
while tracking down the right side of the potential be- 
cause there is no tracking solution to the field equations 
for wbg < 0. 

Once <j) reaches the minimum of the effective potential, 
4>mini it evolves with it. 4> m in increases with expansion 
because Vm is decaying — that is, the interacting dark 
matter is diluted by expansion. (f> can only roll down 
Vq as fast as (j> m in evolves, and for steep potentials it 
acquires a potential dominated state with the following 
properties: 
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where pq = p — — pm is defined as the quintessence 
energy density, and pu = Vm the total matter energy 
density. 

Therefore 4> is in the minimum today, evolving to dom- 
inance with a current equation of state wqq = — 1 + 
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Figure 1: The dashed line represents the <p self-interaction 
potential, the dot-dashed line the effective potential from the 
interaction with iCDM (which decays with expansion), and 
the solid the sum. At early times <j> tracks down either side 
(either A\ or A2), coming to and then slowly rolling with the 
minimum (B). 



3/Qqq (a + f3) . In the future, the Universe is asymp- 
totically approaching a state where baryons and nCDM 
have been diluted away leaving only quintessence and in- 
teracting dark matter in the ratio pu — P%cdm = jjPQ- 
Note that iCDM becoming significant near the onset of 
quintessence domination is not an arranged coincidence, 
but rather an unavoidable feature of this model due to 
the interactions. When <j> is at the potential minimum, 
as it is during quintessence domination, the interaction 
guarantees that the quintessence and iCDM densities are 
comparable. 

One might be tempted to omit the degrees of freedom 
which we label nCDM |l8|. However, without dilution 
by nCDM, the energy pumped into the iCDM causes 
the matter to dominate over many fewer e-folds and less 
time and would drastically alter the matter power spec- 
trum (since matter-radiation equality would be too re- 
cent). Furthermore, the ratio (3 /a would be set by the 
present ratio of energy in CDM to baryons (~ 5 : 1) to 
approximately 6£Iqq/5Qmo- In order to have critical den- 
sity at the present temperature of the cosmic microwave 
background, a and (3 would have to be tuned to values 
that imply a wq today inconsistent with observations. 
Our model just requires an inequality between them to 
be satisfied. The amount of nCDM needed to dilute the 
energy pumped into iCDM depends on the value of (j) c , 
but the amount needed to allow (3 /a to vary is set by 
PnCDM ~ (1 - 6aftQ /5(3n M o) Pcdm- 

For our model to address the coincidence problem, 
achieving a constant ratio of dark energy to dark mat- 
ter densities in the future is not sufficient. It must be 
"natural" for quintessence to dominate only after matter- 
radiation equality, though Vq may be due to physics at 
much higher scale. This yields a mild lower limit on K<p c : 
K(f) c > ^ln(m Q /p eq ). For m Q ~ TeV, k(/) c > 108/a. 
Large a reduces the size of the required <p c and is also 
desirable for a large negative late-time quintessence eqUa- 
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Figure 2: The fraction of the energy density in radiation, mat- 
ter and quintessence as a function of the natural log of the 
scale factor N. The thin dashed line is the tracking attractor, 
while the thin dotted line represents the minimum of the ef- 
fective potential. Note <fi follows the former until it crosses the 
latter. For this numerical example a — 5,/3 = 15, Qqo = 0.6. 
Note that in the future, the ratio of the dark matter and dark 
energy densities approaches the constant a//3. 



tion of state. Thus we see that the same feature of this 
model — large a, (3 with (3 > 2a — simultaneously avoids 
early cosmic acceleration both while <j) is tracking and in 
the minimum of V, as well as yielding better agreement 
with <1l(z) inferred from SNIa data, as we shall show 
now. 

For connecting to observations it is useful to be able 
to talk about the quintessence and the iCDM compo- 
nent separately. Due to the coupling, it is somewhat 
arbitrary how to split the <f> energy density into the 
quintessence part and the interacting dark matter part. 
The choice will lead to correspondingly different pres- 
sures in these components. Of course, for predictions of 
actual observables the choice does not have any effect. 
Thus, for comparison with non-interacting quintessence- 
CDM models, we use the split implied by the above ex- 
pression for Wtot- With this convention, the equation 
of state of quintessence today can be arbitrarily close 
to that of a cosmological constant. Alternatively, if in- 
stead of pressure one uses the decay rates as a guide for 
how to split quintessence and iCDM and assign equa- 
tions of state, then both quintessence and iCDM have 
w = —[31 (a + f3), while the non-interacting matter has 
W = 0. Since SNIa data probes cIl(z), an integral of the 
scale factor, it 'sees' the iCDM as if it were part of the 
quintessence, both decaying as if w = —(3/ (a + (3). Mea- 
surements of the matter density based on gravitational 
dynamics 'see' piCDM a s matter. Thus SNIa and large- 
scale structure measurements would see a different £Imo 
in this model, a possible explanation of their disagree- 
ment for wq > — 1 which was noted in pj. 

A surprising feature of iQCDM models is that they can 
predict cIl(z) that are uniformly larger than predicted in 
AC DAI models, "simulating" w < — 1 at the current level 
of the SNIa data. The energy transfer from quintessence 
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Figure 3: Luminosity distance-redshift (cIl(z)) curves for 
ACDM, a "Big Rip" equation of state wq = — 1.4 (allowed by 
current SNIa data 0) and iQCDM (a = 5,/3 = 15,«qo = 0.6). 
The inset shows the percentage difference from the ACDM 
curve. Also shown in the inset (dot-dashed) is the curve one 
would obtain if the quintessence and iCDM were both treated 
as fluids with equations of state given by their decay rate: 
w = -/?/ (« + /?) = 0.75. 

to a species of matter causes that matter to have a slower 
decay rate. Thus one expects a larger luminosity distance 
for a given redshift than one would obtain from a model 
where quintessence did not interact and had a constant 

As shown in Fig.|3for an example with (a = 5, /3 = 15, 
^qo = 0.6), the current equation of state of quintessence 
is close to —1, and thus <1l{z) is similar to a model with 
the same JImo and wq less than —1 (though the shape is 
different), or nearly identical to a non-interacting model 
with n Q0 = Qqo (0 + a) //3, and w Q = -0/ (/3 + a). As 
SNIa observations improve in both number and redshift 
depth, if we measure a larger luminosity distance than in 
a ACDM model we have the choice between two possible 
conclusions: pq increases with expansion as for example 
in theoretically problematic phantom energy models or 
quintessence and dark matter interact. If we can deter- 
mine Qmo by other means we may be able to rule out 
one or the other based on the shape of (Il(z). Since the 
physical equation of state never drops below w = — 1, 
objects that are gravitationally bound today will remain 



so forever. 

A very important general aspect of iQCDM models is 
that the iCDM energy density changes differently with 
expansion to pure dilution. In this example, when <p 
is in the potential minimum, energy is transfered from 
quintessence to the iCDM. As the quintessence comes 
to dominate at late times, the iCDM becomes signifi- 
cant, and pm {scale factor) 3 increases. This could cre- 
ate observational signatures which have not been looked 
for thus far. It would thus be useful to reexamine pre- 
vious experimental determinations of the matter den- 
sity to place constraints on pm (scale factor) 3 at dif- 
ferent epochs. Large-scale structure data must be re- 
considered: because pu [scale factor) 3 is not constant, 
one must carefully consider at what epoch an observa- 
tion measures it. Energy transfer and the fifth-force will 
have significant effects on structure formation in iQCDM 
effects. Detailed calculations of structure formation in 
these models are ongoing. Preliminary results show that 
there are combinations of the parameters in our model 
that fit observations of the cosmic microwave background 
as well or better than standard models pH ]. 

To summarize: we present a class of models where we 
allow couplings of the dark energy with dark matter de- 
grees of freedoms. Within these models the same con- 
dition that avoids early cosmic acceleration also makes 
the quintessence equation of state more negative today. 
Note that the parameters of the Lagrangian only need 
to satisfy inequalities for these results to hold, obviating 
the need for fine tuning. If our dark sector is iQCDM- 
like, we would predict a larger luminosity distance to a 
given redshift than for non-interacting dark energy with 
the same constant equation of state. Therefore current 
indications of w < —1 can be re-interpreted in our sce- 
nario, removing the impending doom of a "Big Rip" in 
the near cosmic future. 
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